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Abstract 

. We present a method of generation of exact and explicit forms of one-sided, heavy-tailed Levy 

stable probability distributions ga{x), 0<x<oo, 0<Q!<1. We demonstrate that the knowledge 
of one such a distribution Qaix) suffices to obtain exactly gap{x), p = 2,3, ... . Similarly, from 
known ga{x) and g/s^x), < a,/3 < 1, we obtain gapix). The method is based on the construction 
of the integral operator, called Levy transform, which implements the above operations. For a 
rational, a = l/k with I < k, we reproduce in this manner many of the recently obtained exact 
results for gi/kix). This approach can be also recast as an application of the Efros theorem for 
generalized Laplace convolutions. It relies solely on efficient definite integration. 
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I. INTRODUCTION 



The one-dimensional Levy stable one-sided distributions ga{x), < a < 1 are normalized 
probability density functions (PDF) defined on < x < oo, with the peculiar property 
that their mean, variance and all higher moments diverge. The index a enters the principal 
characterization of ga{x) via their Laplace transform 

POO 

^[9a{x);p\= e-P^ gaix) dx = e-P" , p > 0, (1) 

^0 

n 

which is sometimes referred to as the Kohlrausch- Williams- Watts function [1]. The diver- 
gence of all the moments of ga{x) of positive integer order is related to a slow decay of ga{x), 
(gaix) — )■ x~^^^"\ as X — )■ oo), usually termed as a "heavy-tail". Theoretical properties o" 
gaixVs are discussed in j2|, whereas the review of physical applications can be found in js] 
and [41. Many applications of 5'q(x)'s in financial mathematics can be found in jsl, whereas 
their use in various simulations of jump processes is exposed^in The quest for exact 



and explicit forms of ga (x) was accomplished only recently [7|-|9| . The universal formulas of 
Eqs. (2), (3) and (4) in 7^ give explicit forms of gaix) for fractional a = l/k with / and k 
arbitrary relatively prime integers, fulfilling the condition I < k. These results are satisfying 
as in the general case they furnish gi/k{x) as a sum of k — 1 hypergeometric functions, which 



is a neat confirmation of the conjecture of Scher and MontroU [10|]. As indicated in p, Isl the 
implementation of Meijer's G functions as well as of the hypergeometric functions in recent 
version of computer algebra systems pJi] greatly facilitates the practical use of gi/k{x) in 
applications. 

The purpose of this work is to present an alternative approach to generate gi/k{x) by 
an iterative procedure based on the self- reproducing property of ga{x), directly resulting 
from Eq. ([1]). We shall proceed in an increasing order of complexity treating first in Sect. 2 
the low-order values of a, (a = 1/2,1/3 and 2/3), for which exact forms were obtained 

n n n n 

in [7|, Il3|, [ij], see also [15|. In Sect. 2 we also formulate the formalism for general a and 
develop the notion of Levy transform (initiated in [l3]), which will allow one the generation 
of general distributions gi/k{x) by simple integration. In Sect. 3 we present applications of 
our method. Sect. 4 is devoted to discussion and conclusions. The definition and notations 
of Meijer's G functions are given in the Appendix A. A compact formulation of our results 
using the Efros theorem is presented in Appendix B. 



II. LOW-ORDER a = l/k AND REPRODUCING PROPERTIES OF 5„(x) 



a) The content of this subsection a) concerns the first genuine stable distribution gi{x) 
discovered by Levy |l2j, and therefore all the results enumerated here are known and are 
included only for sake of completeness. This historically first Levy stable distribution is: 



3 ; 
2 



X > 0, 



(2) 



satisfying £ 



e with p > 0. We can formally introduce an additional variable 



t in the problem by writing 



X 



-iV(4x) 



(3) 



2^/^T X2 

This last function is of particular interest for us as it intervenes directly in the follow- 
ing property of standard Laplace transform (see formula 1.1.1.26, p. 4, vol. 5 of tl6i]): if 
£ [f{x);p] = F{p), then if we define such that 



^ [ t exp (-^ ] f{t)dt, 

2J^xl Jo V 4x,' 



then 



hix);p =F{p^). 



The proof of Eq. follows: 







POO 


1 r 

or'- 


jC 


h-{x);p 


= / e-P" 




Jo 





f{t)dt 



dx 



(4) 



(5) 



(6) 



2v^ 



fit)t 



e e 



X2 



dt. 



(7) 



To evaluate the internal integration in Eq. ([7]) we use the formula 2.2.2.5, p. 31, vol. 4 of 
JigI . and Eq. ([7]) becomes 



and therefore, from Eq. 



f{t)dt = F{p-2), 



(8) 



£ 



t exp ( 



2y/nx 



3 ' 
2 



p 



-tVP 



(9) 
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We note that although all our PDF's are one-dimensional, gi{x/t^) is formally equal to 

the three-dimensional heat kernel, through appropriate renaming of the variables. 

b) As a next s tep in our development we consider the stable Levy density for a = 1/3 whose 



exact form is 
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15| 



1 



X > 0, 



(10) 



37rx2 ^ V3v3x, 

where Ki,{x) is the modified Bessel function of the second kind, which can also be expressed 
through the Airy function Ai{y) giving 

Ai f (3a;) 



(3x^ 



(11) 



We shall evaluate the Laplace transform of gi{x) using the formula 3.16.3.7, p. 355, vol. 4 
of [l^, that we quote here |l7|: 



X 2 Ki 

3 



27r 



exp 



2\ 3 



4 



p>0, 



(12) 



which, upon choosing 3(a^/4)3 = t, transforms into 



-(-)'''<■ 

37r \xj 3 \ 3^3; 



■,P 



-IP'S 



or equivalently, can be rewritten as 



X 



(13) 



(14) 



In full analogy with Eqs. 
/i(x) defined through 



and ^ we conclude that if £[/(x);p] = F{p) then for 



1 



37rx2 



3 / 2ti 
3 \ 3V3x 



C[fi{x)]p] = F{p-3 



f{t)dt, 



x>0, 



(15) 



(16) 



c) We continue our review of Levy stable densities expressible with standard special functions 
with g2{x) which reads 

2^3 



g2[x) 



277rx3 



2 

e 27 0:2 



Ki 



3 \27x' 



+ K2 



3 \27x^ 



x> 



(17) 



X 3 



5/6 
2/3 



S'^x"^ 



2/9 _r / 7/6 



33 



X 



2 ' 



(18) 
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where in Eq. f|T8|) = iFi{a;b;z) is Kummer's confluent hypergeometric function. 

In analogy witfi Eqs. ([3]) and f|T4|) we form t~2 (^2 (x/tz) in the version of Eq. (fT8|) : 



1 /x 



£(1)4^^/5/6 

a/Stt \ 2/3 



2+3 



^^^^^^ / 7/6 
r (I) a;i ^ M 4/3 



-2H 



2+3 



33x2 



(19) 



To calculate the Laplace transform of Eq. (1191) we apply the formula 3.35.1.16, p. 511, vol. 
4 of jl6| and, after simplifications, we arrive at the identity 













(i) 




= e 











This yields in turn, for £[/(x);p] = -F(p), the pair of equations 

1 fx 



(20) 



and 



/|(x)= / -giA-] fit)dt. x>0 

\t2 



C[f2{x);p] = F{pr 



(21) 



(22) 



in complete analogy to the pairs of Eqs. (jlj), and Eqs. (IT31) . (fTB]l . 

d) The above pattern suggests the validity of the following general integral identities for 
£[/(x);p] = F{py. if 



\9a(^] f{t)dt, 0<a<l, x>0, 

ta \ta 



(23) 



then 



£[/„(x);p] = F(p"). (24) 
We are going to prove Eq. using Eq. ([1]). To this end we calculate explicitly Eq. I^^: 



-px 



^da ) f{t)dt 

ta \t 



tc 



m 
fit) 

POO 

/ f{t)e-P'''dt = F{p'') . 
Jo 



dx 



dt 



dt 



(25) 
(26) 
(27) 
(28) 



In Eq. (126|) we have used a simple change of variable, whereas in Eq. (1271) we have applied 
Eq. (dD. 

Eqs. (123 p and f l24p acquire additional importance when for f{t) we choose another Levy 
stable density, say Qisix) with arbitrary /3, such that < /3 < 1. Then 

/ ( 4 ) 9p{t) dt= [ ( 4 ) g^{t) dt, a; > 0, (29) 

Jo t- J Jo ti^ J 

which defines the following transitive property of Levy laws: 

9ap{x)= I \9a\^] 9ii{t)dt= f \9p\^] 9a{t)dt, (30) 
Jo tc. Vt"/ -'0 t/3 Vt/s/ 

for < < 1. Eqs. and (1301) constitute the key results of the present investigation. 

The following remarks are in order. 



i) The formulas Eqs. (jlj) and ([5]) deserve to be better known 
their several extensions) in various works of Russian school 
spicuously absent from other monographs and handbooks 



he y a ppear (along with 



16 



19 



22 



20|], but are con- 



ii) the formulas Eqs. f|T6|) . fl20|) and more generally, Eqs. fl23l) and fl24l) should complement 
the formula of Eq. ([5]) in the lists of properties of Laplace transform. In subsequent 
Sections we obtain many other formulas of this type, by specifying other values of a 
and (3. 

iii) Eq. (!30|) should be viewed as a tool to generate the Levy stable PDF's, starting with 
those with arbitrary a and /3 and yielding one with a /3, via integration. 

The content of Eqs. fl23|) and fl24l) can be seen as an integral transform with the positive 
kernel Ka{t, x) = t~ag^{^x/ta) with x, t > 0, which we define as the Levy2 transform of index 
a, and is denoted by La'': 

POO _ 

L'i^[fity,x]= Ka{t,x)f{t)dt=~Ux). (31) 

^0 

According to Eq. ( I30l) the action of La on Levy's PDFs satisfy 

L^^\9p{t)-x\ = Lf[ga{t)-,x] = gap{x), (32) 

which well illustrates the reproducing property of Levy densities under the Levy2 transforms. 
It can be seen from Eqs. (|3T1) and (|32|) that for < a < 1 the inverse of Levy2 transform 



cannot be meaningfully defined. As a consequence, Eq. (132|) reflects the semigroup property 

/r)\ 

of La in accordance with the general theory [26j . 



Alternative way to derive the results of this section is to use the Efros theorem of theory 
of Laplace transform, see Appendix B. 



III. SOME APPLICATIONS OF THE METHOD 



a) We shall exemplify now how the method works by choosing in Eq. fl32|) a = /3 = 1/2 
which accordingly should furnish gi{x). We use Eqs. (E]) and ([3]) and obtain 



t 2 exp 



\Ax ^ At 



dt. 



(33) 



47r X2 Jo 

The integral in Eq. (!33|) is not elementary, but with a new variable y = it becomes 



a particular case of the formula 2.3.2.14, p. 322, vol. 1 of [23| (the alternative way is to 
employ the formula 2.2.2.7, p. 32, vol. 4 of 16|), which involves a finite sum of hyperge- 
ometric functions. The final answer, after the simplifications in hypergeometric functions 



z) = pFg (ai, ...,ap]f3i,..., (3g; z), is: 



r(!) 



- 0-^2 



5 3 
4' 2 




4i/vr X 



- 0-^2 



3 5 
4' 4 




+ 



4r(!) 



2' 4 



(34) 



which is precisely the result obtained in j3] a nd ^|13| . We stress that here it was calculated 
with one definite integration using the tables [23]. 

b) Next example in this spirit will involve the indices a = 1/2 and /3 = 1/3, i. e. we shall 
explicitly perform the Levy2 transforms of type 



Lf[gUty,x] = L'^>[gdty,x] = grjx), 



-(2)r 



with Eqs. ([3]) and f|TO|) . The result reads (6 = 

, . 1 



6 



2y/7l X 



3 

2 Jo 



te 



(35) 



(36) 



which, with a change of variable, transforms into 



1 



6 



12('7rx)2 







y ■* Ki 



dy. 



(37) 



The integral in Eq. fl37|l can be perceived as the Laplace transform of y~iKi{by~^^^) for 
which an exact formula in terms of a specific Meijer's G function exists, compare Eq. 3.16.3.9, 
p. 355, vol. 4 of |l6|, with the specializations: a = b = I = 1, k = 2, fi = —3/4, 
u = 1/3 and p = l/(4x). The final result is proportional to the Meijer's G function of type 

0,5 



G 



5,0 



■ X 



A JL 3 11 13 
12' 12' 4' 12' 12 



see Appendix. By lumping together all the constants we obtain: 



gi{x) 



V2 



A X 3 11 13 
12' 12' 4' 12' 12 



(38) 



¥m 



2-^0-^4 

X6 



112 5 
3' 2' 3' 6 




6r(i) 



— 0-^4: 



X3 



12 5 7 
2' 3' 6' 6 



H 0-^4 



TT X2 



2 5 7 4 
3' 6' 6' 3 




V3r(|) 

727rx3 



o-fk 



5 7 4 3 
6' 6' 3' 2 




3-i[r(i)]' 

^ „ 12 3 11 0-^4 

2 3 7f2X a 



7 4 3 5 
6' 3' 2' 3 




(39) 



The Eq. (!39|) was obtained using the formulas 16.17.2 and 16.17.3 of [24]. The result of 
Eq. ( l39|) is a special case of the solutions in The exphcit form of gi{x) in Eq. ( l39|l 

6 

appears to be written down for the first time here. 

c) The third example in this section will be of more academic value as it will concern three 
PDF's known from the previous sections: gi{x), gi{x) and 5'i(x). According to Eqs. ( I30|l 
and fl32|) they are related through 



Lf[g.{t)-x] = Lf[gAt)-x]=gAx). 



(40) 



8 



Indeed, the first integral transform of Eq. (HO!) is equal to 



45 3 



t2 



r(i) 



2^3 



TT X3 Jo 



t 2 e 4a; 




2+3 



9v^r (I) x7/3 



t2 e 4^ 



and, with the help of the formula 3.38.1.30, p. 553, vol. 4 of [16j, it can be written down as 




dt. 



(41) 



a special case of the Mel 
to Eq. fl38l). however see 



er's G function, very much in the spirit of considerations leading 



25[. (Alternatively the simplified formula 3.35.1.16, p. 511, vol. 4 



of [l6| can also be used.). We will skip further details of this evaluation and conclude that 
the result is equal to gi{x) from Eq. ffTOj) . as it should be. 

d) In the following examples we shall not use the reproducing property of Eq. (l30ll . We 
shall instead concentrate on Eqs. fl2^ and fl2^ conceived as an operational property of the 
Laplace transform. To this end we, quite arbitrarily, choose the modified Bessel function 
Ko{x) as /(x) in Eq. (!23l) and we read off its Laplace transform from the Eq. 3.16.1.2, p. 



349, vol. 4 of 



16|, which is 



c[Mxy,p] 



arccos(p) 



1 _ p2 



p> 0. 



(42) 



If we set gi {z) in Eq. ( I23ll then the corresponding transformed function is 

e^r(0,a;) 



^0,1 (^) 



TTX 



e^Ei{—x) 
2^/ttx 



(43) 



(44) 



In Eq. ( H3ll we have used the formula 2.16.8.5, p. 352, vol. 2 of |23[, with r(0,x) = — 
Ei{—x), see p. 726, vol. 2 of [23|. In Eqs. ( H3|) and f H4|) . r(z/, z) and Ei{z) are incomplete 
gamma function and the exponential integral, respectively. The Laplace transform of the 
function in Eq. (144|) can be evaluated to be 

arcsinh ( ^/p^^] 



Knl{x)]P 



Let us briefly sketch how Eq. (H5|l comes about. 



p > 0. 



(45) 



We write out explicitly the Laplace transform in Eq. (jl]) as 



2\/ttx 



dx 



' 1 p-1 



2' 



P 



arcsinli (^y/p — l) 

Vp^ 



, p>0. 



(46) 



(47) 



(48) 



In Eq. f H7|) we have used the 
Eq. 7.3.2.83, p. 473, vol. 3 of 



in. 3.4.1.3, p. 135, vol. 4 of [16], whereas Eq. fHS]) results from 
23| . Since for p > 0, arccos(^) j \/\ — p = arcsinh(A/p^n^)/A/j)^n^, 



the above calculations evidently confirm the validity of Eq. (|2i 
We go now a step further: while still keeping f{x) = Kq{x) in Eq. (^^, we use now gi{z), 
see Eq. (fTO!) , to perform the integration in Eq. (!23l) . The corresponding transformed function 
Kq 1 (x) reads |ll | 



^3/2 



(-J /I R I ^ 



4^3/2 4/ 



) TC^X' 



_5_ X ii 13 

12' 12' 12' 12 

A A 3 3 13 13 
12' 12' 4' 4' 12' 12 



X > 0, 



(49) 



whose Laplace transform can be calculated via 3.40.1.1, vol. 4 of 16|. This lengthy expres- 
sion involves a combination of four different hypergeometric functions of argument p^. We 
shall not reproduce it here. We emphasize however that we have verified numerically the 
relation 

arccos (p-*^/^) 



Koiix);p 



1 — p 



,2/3 



p> 0, 



(50) 



thereby establishing an algebraic identity stating that the aforementioned combination of 
hypergeometric functions is equal to the r.h.s. of Eq. (I5U]) . 

It is clear that these procedures can be carried out by choosing functions other than Ko{x) 
in Eq. ( |23l) . This circumstance paves the way for a scheme of generation of algebraic iden- 
tities involving special functions using Eqs. (!23|) and (^^. This will be a subject of future 
publication. 
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IV. DISCUSSION AND CONCLUSIONS 



Our kernel Ka{t, x) which defines the Levy2 transform of index a through Eq. (13T]) above, 
can be related to the function n{s, r) introduced by E. Barkai in jisj] as follows 



r9a[ — ] = — fia{r,s), 

a s \sq / as 



(51) 



where attention should be given to the order of variables in Eq. f lSTj) . This function serves as 
a kernel in another integral transform 13| which links a solution of ordinary Fokker-Planck 
equation Pi(x, r), satisfying, for given functions A{x) and B{x), 



-— [Aix)P^ix,r)] + — [B (x) P,ix , r)] 



(52) 



with Po(x,r) satisfying the fractional Fokker-Planck (FFP) equation 



^ [A{x)P^{x, r)] + ^ [Bix)P^ix, r)] 



(53) 



where is appropriately defined fractional derivative. The relation is 



Paix,T) 



n{s, t) Pi{x, s) ds, 



(54) 



and is called in 



13| the inverse Levy transform. The functional relation of this type appears 



also in Sokolov's work [27|, section III, Eqs. (15) and (16) of [27]. If we provisionally reserve 



the notation ( L 



'(1) 



to the transform in Eq. ( l54l) . then it can be rewritten as 
{L^^^y\p^{x,s)■r]=P^{x,T). 



(55) 



The symbol La^ should be clearly distinguished from the Levy2 transform appearing in 
Eqs. ( pTl) and ( l35ll . denoted by Lq, . From Eq. (|5Ti) the relation between those two transforms 
is 



1 



a s 



-Pi(x, s);r 



(L«)-^[A(x,s),r] 



(56) 



where it should be stressed that in both sides of Eq. (1561) the integration variable is s. Here 
for the general a the direct Levy transform La^ cannot be constructed, compare remarks 
after Eq. (ES]). 



In both Refs. 13|] and [27| the emphasis is on generation of solutions of FFP equations 
depending on the choice of A{x), the drift function, and B{x), the diffusion function, whereas 
our method here focuses on generation of explicit forms of ga{x). In the separate work we 
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have undertaken the analysis of FFP equation using the exphcit form of gi{x), see [281. 
Further considerations on apphcations of relations of type Eq. f l5^ can be found in 29 1 
and [30]. Extensions of the current framework to include two-sided Levy distributions 8|, |9| 
as well as log-Levy distributions 3l| are under study. Finally, note that newer versions of 
Mathematica® furnish numerical forms of ga{x)- Similarly, Matlab® can be also used for 
that purpose. 
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Appendix A: Definition and notations of Meijer's G functions 



The Meijer's G function is defined as an inverse Mellin transform denoted by ^, see 
vol. 3 of ^]: 



G 



PA 



(ap) 



M 



-1 



nr=ir(/3.+^)n;=ir(i 



aj — s) 



n 



j=m+l 



r(i-/3,-^)m=.+irK + .) 



; z 



G{[[ai, ...,an], K+i, • • • , ftp]], • • • , /3m], [Am+l, • • • , f3q]], z), 



(Al) 



(A2) 



where in Eq. ( lAip empty products are taken to be equal to one. In Eqs. (1A1|) and (1A2|) the 
parameters are subject of conditions: 

z 0, < m < q, < n < p; 

aj e C, J = 1, . . . ,p; /3j G C, J = 1, . . . , g; 



(ap) = ai, a2,..., ap] {(3g) = f3i, (32, ■ ■ ■ , (3q. 



(A3) 
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For a full description of integration contours in Eq. ( lAip . general properties and special 
cases of the G functions, see vol. 3 of [23] • In Eq. flA2|) we present a transparent notation 
inspired by computer algebra systems 

In this context we can give more details on how Eq. fl5S]) comes about. The rewriting 
again of Eq. 3.16.3.9, p. 355, vol. 4 of [16] with the set of parameters given above, gives the 
following expresion for (71 (x): 



G 



0,5 
5,0 



_5_ X 3 11 13 
12' 12' 4' 12' 12 



(A4) 



According to (IA1|) . in Gjq the parameters are: m = 0, n = 5, p = 5, and g = 0. If we apply 
this to the notation specified in the function Gg'g in flA4l) . it becomes: 



G 



0,5 
5,0 



6^x 



_5_ X 3 11 13 
12' 12' 4' 12' 12 



G 



5 7 3 11 13 
12' 12' 4' 12' 12 



which is a required format of computer algebra systems 



[[]'[]]'6^^ 



(A5) 



All the Meijer's G functions encountered in course of integrations via Levy2 transform of 
Eq. fl3T]) were converted to forms of type flA5|) and consequently checked numerically. 



Appendix B: The Efros Theorem 



We quote here the A. M. Efros theorem (1935) which is the generalisation of the con 



19 



32 



volution theorem for the Laplace transform. See Refs. 
applications. 

Theorem: If G{p) and q{p) are analytic functions, C [/(x);p] = F{p) and 

/"OO 

C[g{x,t)-p\ = / g{x,t)e~^-dx = G{p)e-"i^^\ 
Jo 



33| for proof and many 



then 



G{p)Fiqip)) 



f{t)g{x,t)dt 



e-P'^dx. 



(Bl) 



(B2) 



We shall demonstrate now that Eq. ( 130|) can be recast as a special case of the Efros theorem, 
Eq. ([B2]). 
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The use of the result (lB2p in the context of Levy stable laws entails the choice f{t) = gp{t), 
< /3 < 1, implying F{p) = e~^'^ with Eq. ([1]). Furthermore, we set G{p) = 1, q{p) = 
and 

1 / X \ 

< a < 1. 



9i^^t) = ^9a (^) 

From the Efros theorem (lB2p we obtain 



G(p)F(g(p))=e-^"' 



On the other hand, Eq. ([T]) implies 



9isit):^9a (^) dt 



e-P'^dx. 



which immediately yields 



9ais{x)= / 9is{t)-^ga[:^) 



dt. 



(B3) 



(B4) 



(B5) 



(B6) 



which is our Eq. (!30!) . 
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